Electron-correlation effects in appearance-potential spectra of Ni 
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Spin-resolved and temperature-dependent appearance-potential spectra of ferromagnetic Nickel 
are measured and analyzed theoretically. The Lander self-convolution model which relates the line 
shape to the unoccupied part of the local density of states turns out to be insufficient. Electron corre- 
lations and orbitally resolved transition-matrix elements are shown to be essential for a quantitative 
agreement between experiment and theory. 



I. INTRODUCTION 

Due to the Coulomb interaction, an electron approach- 
ing a metal surface can excite a core electron into a 
state of the unoccupied bands. Above a threshold en- 
ergy both, the excited core electron and the de-excited 
primary electron are scattered into valence states just 
above the Fernii, ejuergy. The appearance-potential spec- 
troscopy (APS)ldoBa monitors the intensity of this radi- 
ationless transition as a function of the energy of the pri- 
mary electrons by detecting the emitted X-rays or Auger 
electrons of the subsequent core-hole decay. Since core- 
hole formation is involved, the method is element specific 
and local. Its comparatively simple experimental setup 
and its surface sensitivity qualifies APS as a useful tech- 
nique for surface analysis. For a ferromagnetic material, 
the spin dependence of the AP signal obtained by using 
a polarized electron beam reflects the surface magnetiza- 
tion as has,|t)|eeit,demonstrated for the transition metals 
Fe and Ni.ErBiDla'Ll In a sense APS can be considered to be 
complementary to Auger-electron spectroscopy (AES).t3 
While the Auger line shape resulting from core-valence- 
valence (CVV) transitions yields information on the oc- 
cupied part of the valence band, APS is sensitive to the 
unoccupied electronic structure. Contrary to k resolved 
(inverse) photoemission, the AP transition is more or less 
localized in real space. This suggests to employ APS as a 
quantitative probe for the unoccupied part of the (spin- 
dependent) local valence density of states (DOS). 

In most cases the interpretation of appearance- 
potential spectra is still based on a simple independent- 
electron model suggested by Lander in 1953:Eil Hereafter, 
the line shape is given by the self-convolution of the 
unoccupied part of the DOS. The desired informatioa. 
can then be obtained by de- convolution techniques £3't3 
Within the context of AES, however, it soon became 
clear that the self-convolution model seems to be over- 
simplified. Powelllia discovered an "anomalous" shape of 



the CVV Auger line of Ag which has been attributed to 
electron-correlation effects. Correlations may be signif- 
icant also for APS because of the direct interaction of 
the two additional valence electrons in the final stata^ 
This is demonstrated by the Cini-Sawatzky theory:t£lE£l 
Within the framework of the single-band Hubbard model, 
the two-particle (APS/ AES) excitation spectrum is dom- 
inated by a pronounced satellite feature which is split off 
the band-like part by a characteristic energy of the order 
of the on-site Coulomb interaction U . For more realistic 
(multi-band) Hubbard models the direct rqrrdaticma give 
rise to very complex satellite structures. OE3ll3c3'El 

A second shortcoming of the Lander model consists in 
the fact that transition-matrix elements are not taken 
into account. Modern theories for APS (AES) based on 
density-functional theory and the local-densityr-aporox- 
imation (LDA) have overcome this deficiency.c3E3 Al- 
though the direct correlations are still neglected in these 
approaches, it seems that the AE|Sipectra of Fe and Ni, for 
example, are well understoodHo - serious indications 
for correlation effects have not been observed. This is sur- 
prising since it is well known that the Auger line shape of 
Ni cannot be exolained within an effective independent- 
electron modelHj It is dominated by a strong spectral- 
weight transfer due to an on-site interaction of the order 
of C/ w 2 — 3 eV tpibe. compared with the effective c?-band 
width A « 3 cV.BeI 

Clearly, the question what really determines the 
appearance-potential line shape of transition metals, can 
only be answered a posteriori - namely by compari- 
son with a theory that realistically includes electron- 
correlation as well as matrix-element effects, orbital de- 
generacy and sp-d hybridization from the very beginning. 
The analysis of the AP line shape of Ni, as a prototypi- 
cal ferromagnetic 3d metal, is the purpose of the present 
paper. In the experiment spin-resolved spectra from ex- 
citation of the 2p3/2 {Liii) core level are recorded for 
different temperatures ranging from T = 100 K up to 
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and slightly above the Ni Curie point Tq = 624 K.EZI 
The theoretical interpretation of the line shape is essen- 
tially based on a multi-band Hubbard-type model with 
general on-site Coulomb interaction and a realistic pa- 
rameterization of the (LDA) one-particle electronic struc- 
ture as an input for a subsequent many-body calcula- 
tion. Standard diagrammatic techniques are used to ac- 
count for correlations. The transition-matrix elements 
are calculated within the usual intra-atomic approxima- 
tion. The novel feature of the approach is that it allows to 
study separately the DOS effect (bare self-convolution), 
the temperature-dependent effect of the (direct and re- 
maining indirect) correlations as well as the effect of the 
transition-matrix elements. 



II. EXPERIMENT 

Experiments are performed for a Nij^iip) single- 
crystal surface with in-plane magnetization.BcS The spin- 
polarized electron beam used for excitation is emitted 
from a GaAs photocathode irradiated by circularly po- 
larized light. The longitudinal spin polarization of the 
emitted photoelectrons is changed into transversal by 
90° electrostatic deflection. To correct for the incom- 
plete polarization of the electrons {P ~ 30%), all data 
have been rescaled to a 100% hypothetical beam polar- 
ization and alignment between the electron-polarization 
and the sample- magnetization vector. In the present 
setup the core-hole decay is detected via soft-X-ray emis- 
sion (SXAPS). The detector arrangement consists of a 
multichannel plate with filters and a Csl layer acting as 
photon-to-electron converter. The APS signal is mea- 
sured as a function of the primary energy of the elec- 
trons. Modulation of the sample potential together with 
lock-in technique is employed to separate the signal from 
the otherwise overwhelming background. For the poten- 
tial modulation a peak-to-peak voltage of 2V has been 
chosen. This value ensures high APS signals for the 
2p3/2 core level without modulation-induced broadening 
effects. Further details of the experimental setup and the 
sample preparation have been described elsewhere. Q 

Fig. d shows the differential AP intensity as a func- 
tion of the energy of primary electrons with polarization 
parallel (minority, J.) or antiparallel (majority, t) to the 
target magnetization. The displayed energy range covers 
the emission from the Lm transition. The Ln emission 
would be visible at higher energies shifted by the 2p spin- 
orbit splitting of about 17.2 eV. 

For T /Tq — 0.16 the system is close to ferromagnetic 
saturation. The appearance-potential spectrum shows 
a strongly asymmetric intensity ratio as well as a spin 
splitting of the main peak at E — 852.3 eV (dotted line). 
Since Ni is a strong ferromagnet, there are only few un- 
occupied d states available in the majority spin channel, 
and thus /| < /| holds for the (non-differential) inten- 
sities. This is the dominant spin effect. The intensity 
asymmetry in the main peak gradually diminishes with 
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FIG. 1: Spin-resolved Lm Ni appearance-potential spectrum 
for different temperatures T/Tc- Data points: measured dif- 
ferential intensity dl /dE as a function of the primary energy. 
For better comparison with theory the same data are shown 
twice (left and right panel) . Lines, left: theory with direct and 
indirect correlations included (ladder approximation). Lines, 
right: indirect correlations included only (self-convolution). 



increasing temperature and vanishes at Tq. Note that 
only the position of the majority peak shifts with T while 
the minority peak position remains unchanged. 

The main peak is related to the high density of states 
at the Fermi energy and is mainly due to d-d character of 
the two final-state electrons. Additional small s-d contri- 
butions are present in the secondary peak &i E k, 859 eV 
as has been concluded from the analysis of the transition- 
matrix elements.BcJ This structure has been identified as 
originating from a DOS discontinuity |deriving from the 
Lj critical point in the Brillouin zone.cl No temperature 
dependence and spin asymmetry is detectable here. 



III. THEORY 

The theoretical approach is based on the usual two- 
step description assupaing the lifetime of the core hole to 
be sufficiently long.EJo The AP line shape is then unaf- 
fected by the cross section for the core-hole decay step 
but solely determined by the excitation step. According 
to Fermi's golden rule, the intensity can be written as: 

/.,..(k||,i?) « Im ^ M^^^;(k||,£;)x 
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For an incoming electron with spin ai =|, J., one has to 
consider the "singlet" transition, i. e. excitation of a core 
electron with spin dc = — cr^, as well as the "triplet" tran- 
sition with Uc — <Ji- Since the spin state of the final core 
hole is not detected, the intensities have to be summed 
incoherently: 



(k|| , E) = h^^, (k|| , E) + (k|| , E) 



(2) 



Above the Curie temperature 1^^ — I-a-i ■ 

The "raw spectrum" resulting from an intra-atomic 
transition at the lattice site i is given bv-the imaginary 
part of the retarded two-particle (ZubarevEj) Green func- 
tion in Eq. (|l|). It describes the correlated propagation 
of the two additional final-state electrons at energy E: 

(3) 

Here (• • •) is the thermodynamical average, [..., ...]_ the 
commutator, and 0{t) = exp{iHt)0 exp{—iHt) is the 
Heisenberg time dependence of an operator O. The 
Green function is expressed in terms of annihilators (cre- 
ators) ci^j^^ which refer to a tight-binding one-particle 
basis. Due to translational symmetry of the fee lattice 
the i dependence is only formal. L — {I, m} charac- 
terizes the angular momentum of the localized Sd, 4s, 
and 4p basis functions. In total there are 2x9 orbitals 



\iLa) 



vac.) per site. 



Because of the definite L character of the basis or- 
bitals, the orbitals on different sites are mutually non- 
orthogonal: {iLa\i'L'a') = S^^FS^^^. There are different 
advantages to develop the many-body formalism for a 
non-trivial overlap matrix S 7^ 1 (see Refs. |3l| , |3^ ): A 
unique decomposition of the spectrum into s, p, and|-6i 
parts is possible, the basis orbitals arc more localized,Eil 
and in the interaction part of the Hamiltonian the on-site 
Coulomb interaction among the 3d electrons can be dealt 
with separately. Furthermore, the same basis states en- 
ter the definition of the transition-matrix elements, and 
their definite angular-momentum character implies help- 
ful selection rules. 

The transition-matrix elements in Eq. (^, 



(k|h^) = 
(i)(2p,a,|(2)(k||i?a. 



ri -r2 



(4) 



are calculated bai_ assuming the transition to be intra- 
atomic as usual£3o'E3 The different wave functions as 
well as the Coulomb operator l/|ri — are expanded 
into spherical harmonics, the angular integrations are 
performed analytically, and the numerical radial integra- 
tions are cut at the Wigner-Seitz radius (see Ref. ^ for 
details). 



Surface effects enter the theory via the high-energy 
scattering state \\inEai). It is calculated as a conven- 
tional LEED stateo with k|| = to describe the nor- 
mally incident electron beam in the experimental setup. 
At kinetic energies of the order of keV, however, the ky 
dependence turns out to be weak. Furthermore, at high 
kinetic energies, multiple-scattering effects are small and 
may be neglected for convenience. 

The (paramagnetic) LDA potential for Ni is deter- 
mined by a self-consistent tight-Jainding linear muffin- 
tin orbitals (LMTO) calculation.ES The 3c?, 4s, and 4p 
valence orbitals \iLa) are taken to be the muffin-tin or- 
bitals (MTO's). The four- fold degenerate 2^3/2 core state 
is obtained from the LDA core potential by solving the 
radial Dirac equation numerically. Its (relativistically) 
large component is decomposed into a (coherent) sum of 
Pauli spinors |2p, (Tc) with CTc =1,1- 

Note that the Green function in Eq. (|l]) generally de- 
pends on four quantum numbers Li-L^. In fact, for the 
present case each term in the sum gives a significant con- 
tribution to la^a-i- The usual characterization of the fi- 
nal state with two quantum numbers [d-d, s-d, etc.) is 
no longer valid if correlations are included. The orbital 
character may change by electron scattering. 

Neglecting electron correlations altogether and assum- 
ing the matrix elements to be constant, Eq. (|^) reduces 
to the Lander model: In this case the intensity I^. is 
given by the following sum of (singlet and triplet) self- 
convolutions: 

I<y, « J2 j dE' pL,aAE')pL,-.,{E - E') 

L1L2 

+ J PLi-'(^')PL2^^iE ~ E') . (5) 

L1L2 

Here pLa{E) = (1 — f[E))pLa[E) is the unoccupied part 
of the L-resolved and spin-dependent density of states 
(DOS) pLa{E) where f {E) is the Fermi function. The 
prime in the second sum in Eq. (|^) excludes the term 
Li — L2, (Tc = cTi, which is forbidden by the Pauli prin- 
ciple. Note that lattice symmetries require the on-site 
{i = i') but off-diagonal (L ^ L') DOS pLL'a{E) = 

(— l/7r)Im((ciLo-; CiL'cr))s to vanish identically. This is 
a consequence of the choice of the (non-orthogonal) basis 
set \iL(j). 

To estimate the significance of correlation effects, the 
two-particle Green function in Eq. (|l|) is calculated for a 
nine-band Hubbard-type model H — Hi^y^x + Efint ~ A_ff 
including 4s, 4p orbitals and correlated 3d orbitals: 

ii' LL'a 
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iaa' nil .. .m4 



(6) 



The hopping parameters T^f,-^ of the one-particle term 



HhuA (and also the overlap 5"^^/^ parameters) are ob- 
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tained from a Slater-Kqster fit to the paramagnetic LDA 
band structure for Ni.Ej Opposed to photoemission spec- 
troscopy, this comparatively simple tight-binding param- 
eterization appears to be sufficient in the case of APS 
since the two-particle spectrum does not crucially depend 
on the details of the one-particle DOS. 

The on-site interaction among the 3d electrons is de- 
scribed by the second term Hmf The Coulomb matrix 
depends on four orbital indices, C/„iim2m3m4, referring to 
the MTO's for 1 = 2. Using atomic symmetries the inter- 
action parameters can essentially be expressed in terms 
of two parameters U and J.E3 Interactions involving de- 
localized s and p states are assumed to be sufficiently 
accounted for within the LDA. For the d states, however, 
there is the well-known double-counting problem: The 
interactions are counted twice, once in i?LDA (i- e. on a 
mean- field level) and once more in iJint- To avoid this 
double counting, a third term AH has been introd uced. 
by which the Hartree-Fock part of i/int is subtractedE3'E£l 

The Hamiltonian H constitutes an involved many- 
body problem. Due to the low density of 3d holes in the 
case of Ni, however, it appears to be-ceasonable to employ 
the so-called ladder approximatioiiE£l which extrapolates 
from the exact (Cini-Sawatzloi) solution for the limit of 
the completely filled band.E3't3 For finite hole densities 
the ladder approximation gives the two-particle Green 
function {{cc;c'c')) as a functional of the one-particle 
Green function ("direct correlations"). The one-particle 
Green function of the type ((c;c^)) corresponds to the 
(inverse) photoemission spectrum which itself is renor- 
malized by i?int ( "indirect correlations" ) . It is calculated 
self-consistently within second-order pcrtuijtiation theory 
(SOFT) around the Hartree-Fock solution. l3 For a mod- 
erate U ^ A and a lowiiple density, a perturbational ap- 
proach can be justified.E3 A resummation of higher-order 
diagrams is-jmportant to describe bound states ( "Ni 6 eV 
satellite" )lZI which, however, are relevant for AES only. 
The numerical values for the direct and exchange inter- 
action, U = 2.47 eV and J = 0.5 eV, are taken from Ref. 
|32| where they have been fitted to the ground-state mag- 
netic properties of Ni. Since spin-wave excitations are 
neglected in the approach, the calculated Curie tempera- 
ture Tc = 1655 K is about a factor 2.6 too high compared 
with the experiment. Using reduced temperatures T/Tq, 
however, the teroeerature trend of the magnetization is 
well reproduced. l3 



IV. RESULTS AND DISCUSSION 

The solid lines in Fig. |l| (left) show the spectra as cal- 
culated from Eq. (|l|) using the ladder approximation. To 
account for apparatus broadening, the results have been 
convoluted with a Gaussian of width a = 0.6 eV (see Ref. 
p3| ). The calculated data are shifted by 852.3 eV such 
that onset of the unbroadened spectrum for T/Tq = 0.16 
coincides with the maximum of Lm emission in the ex- 
periment (dotted line). 



What are the signatures of electron correlations? To 
estimate first the effect of the direct interaction between 
the two additional final-state electrons, Fig. |^ (right) 
also displays the results of the self-convolution model 
for comparison (still including matrix elements as well 
as the fully interacting one- particle DOS). Fig. |l| and 
also a more detailed inspection show that the secondary 
peak at E Ki 859 eV is not affected by correlations at 
all. This is consistent with observed temperature inde- 
pendence of the peak and with the fact that the DOS 
has mainly s-p character at the discontinuity deriving 
from the L7 critical point. The maximum of the sec- 
ondary peak is therefore used as a reference to normalize 
the measured spectra for each T. Looking at the results 
of the ladder approximation, the overall agreement with 
the measurements is rather satisfying. Except for the 
lowest temperature the intensity, the spin splitting and 
the spin asymmetry of the main peak are well reproduced 
and, consistent with the experiment, a negligibly small 
intensity asymmetry and spin splitting is predicted for 
the secondary peak. Switching off the direct correlations 
(Fig. |l|, right) results in a strong overestimation of the 
main peak structure. Within the Cini-Sawatzky theory 
this has a plausible qualitative explanation: For low hole 
density the main effect of the direct correlations is known 
to transfer spectral weight to lower energies inaccessible 
to APS. This weight shows up again in the Auger spec- 
trum. Recall tijAt-in fact a considerable weight transfer 
is seen in AESLj'E£I and that hypothetically for C7 1— > 00 
all weight would be taken by a satellitp igulit off at the 
lower boundary of the Auger spectrum 

The indirect correlations manifest themselves as a 
renormalization of the one-particle DOS. Thus, in first 
place they are responsible for the correct temperature 
dependence of the intensity asymmetry of the main peak 
in the AP spectrum. In particular, the indirect correla- 
tions result in a narrowing of the DOS, which near the 
Fermi energy is given by the (almost spin indepen- 
dent) quasi-particle weight z ~ 0.88 < 1, and in an in- 
trinsic lifetime broadening F cx {E — Ep)'^,T^. For T = 
the latter turns fowt to be smaller than F « 0.01 eV for 
E — E-p < 1 eV.t3 The discrepancy between experiment 
and theory for T/Tq = 0.16 is possibly due to an un- 
derestimation of F within the SOFT: With decreasing T 
the increase of the exchange splitting leads to the appear- 
ance of a strong peak just above E-p in the minority DOS. 
A larger lifetime broadening of this peak would imply a 
less pronounced feature in the differential AP spectrum 
- mainly in the minority channel. 

Setting M^jJ;(k||,£;) = ±1 = const for Li > L2 or 
Li < L2, respectively, (see below) and comparing with 
the results of the full theory, demonstrates the impor- 
tance of the transition-matrix elements. Their energy 
dependence (via \k^^Eai)) is weak over a few eV at ener- 
gies of the order of keV and cannot explain the difference 
between (a) and (b) in Fig. |^. The main difference is 
rather a consequence of the fact that the radial 2p core 
wave function has a stronger overlap with the (more lo- 
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FIG. 2: Ni AP spectrum for T = 0. (a) full theory, (b) as 
(a) but matrix elements taken to be constant (see text). 



calized) 3d as compared to the (more delocalized) 4s/4p 
radial wave functions. This impUes a suppression of the 
s-j) contributions to the orbital sum in Eq. (Q) as al- 
ready found by Ebert and PopescuO The features above 
E — 860 eV originate from additional discontinuities of 
the s-p-like DOS (as for the peak at « 859 eV). 

For T < Tq the spin asymmetry of the spectrum is 
mainly due to the spin dependence of the Green function 
in Eq. (|^). If the calculation of the matrix elements (^) 
starts from the spin-polarized L(S)DA potential, an ad- 
ditional spin asymmetry is observed resulting from the 
spin dependence of the states in Eq. This, however, 
is small and has practically no influence on the results. 
On the other hand. Fig. g shows a strong suppression 
of the intensity asymmetry at high energies when taking 
matrix elements into account. This effect is controlled 
by the symmetry of the matrix M^j^^ = ^lI'l2^^\\' 
In the antisymmetric case, Ml^l^ = —M^^^ {Li ^ L2), 
there is a maximum spin asymmetry (Fig. g) while, even 
for a ferromagnet, there is no spin asymmetry at all for 
the symmetric case. The results of the full calculation 
along Eq. (|^) are neither fully symmetric nor antisym- 
metric with respect to Li, L2. 



plex interplay of different factors. The present study 
has shown that a quantitative theoretical analysis of 
the temperature- and spin-dependent spectrum must be 
based on three ingredients at least: 

(i) "indirect" correlations: A realistic Hubbard-type 
model including s-p like states is a proper starting 
point to describe the magnetism and the temperature- 
dependent renormalization of the one-particle DOS. 

(ii) "direct" correlations: While s-p derived features at 
higher energies appear to be sensitive to the geometrical 
structure only, the main peak is strongly affected by the 
direct interaction between the two additional final-state 
electrons. Consistent with the Cini-Sawatzky model, the 
dominant effect is a considerable spectral-weight transfer 
to energies below the threshold. 

(iii) matrix elements: Spin-resolved APS cannot 
be described theoretically without calculating orbital- 
dependent matrix elements. The spin asymmetry is 
mainly determined by their transformation behavior un- 
der exchange of the orbital quantum numbers Li L2. 

An open question concerns the importance of core-hole 
effects in APS. Future work shall be concerned with the 
scattering at the core-hole potential in the final state and 
shall include the edge effectg23 known from the studies of 
simple metals. A generalized ladder approximation in- 
cluding the valence-core-jpteraction has been proposed in 



Ref. |20l Previous workliaElO has shown that one should 
not expect the effects due the scattering at the core-hole 
potential to be strong for the case of Ni. However, core- 
hole effects will become more important for systems with 
a smaller 3d occupancy. For Co and Fe one also expects 
stronger effects of d-d correlations. This work has shown 
that these d-d correlations cannot be neglected even for 
a system with low d-hole density such as Ni and that de- 
convolution techniques to extract the unoccupied local 
DOS have to be questioned seriously. 



SUMMARY 
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